Electrochemical synthesis and superconducting phase diagram of Cu x Bi 2 Se3 
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The superconducting Cu x Bi2Se3 is an electron-doped topological insulator and is a prime candi- 
date of the topological superconductor which still awaits discovery. The electrochemical intercalation 
technique for synthesizing Cu x Bi2Se3 offers good control of restricting Cu into the van-der-Waals 
gap and yields samples with shielding fractions of up to ~ 50%. We report essential details of this 
synthesis technique and present the established superconducting phase diagram of T c vs x, along 
with a diagram of the shielding fraction vs x. Intriguingly, those diagrams suggest that there is a 
tendency to spontaneously form small islands of optimum superconductor in this material. 

PACS numbers: 74.25. Dw; 74.62.-c; 82.45.Vp; 65.40.gk 



I. INTRODUCTION 



Topological insulators (TIs) are attracting great inter- 
est since they realize a new state of matter; namely, the 
bulk of such materials is insulating, but in contrast to 
conventional band insulators, their bulk wave functions 
exhibit a non-trivial Z2 topology leading to gapless and 
hence conductive surface states i 3 Those topological sur- 
face states are interesting because they exhibit a Dirac- 
like energy dispersion (similar to that in graphene) and 
a helical spin polarization, both of which hold promise 
for various energy-saving device applications^ — Soon af- 
ter the theoretical predictions of candidate materials^ 
a series of promising materials were experimentally dis- 
covered to be three-dimensional (3D) Tls^r— stimulat- 
ing the search for a superconducting analogue, i.e., a 
topological superconductor, which is characterized by a 
full energy gap in the bulk and the existence of gapless 
surface Andreev bound states^— Such a topological 
superconductor is predicted to be a key to realizing a 
fault-tolerant topological quantum computing^ and the 
discovery of a concrete example of the topological super- 
conductor would have a large impact on future technol- 
ogy. In this context, there is an intriguing prediction that 
when superconductivity is achieved by doping a topologi- 
cal insulator, it may realize a topological superconducting 
stated 

Recently, Hor et al 29 reported superconductivity in 
the electron-doped TI material Cu^E^Ses for 0.1 < x < 
0.3, where the layered TI compound Bi2Se3 was interca- 
lated with Cu. Their samples showed superconductivity 
below critical temperatures T c < 3.8 K. In a subsequent 
study on these samples, 30 it was further found that the 
topological surface states remain intact upon Cu inter- 
calation. Thus, Cu ;E Bi2Se3 is the first promising candi- 
date material to be a topological superconductor, and 
it is very important to elucidate the nature of its su- 
perconducting state. However, the samples prepared by 
Hor et al. showed superconducting shielding fractions of 
only less than 20% and the resistivity never really disap- 
peared below T c j 29 i 3Q Therefore, some doubts remained 



about the bulk nature of the superconducting phase in 
Cu :c Bi2Se3, and preparations of higher-quality samples 
were strongly called for. 

In our recent experiments^ we succeeded in synthesiz- 
ing Cu ;E Bi2Se3 samples with shielding fractions of up to 
50% (depending on the Cu concentration) by applying 
a different sample preparation method than the stan- 
dard melt-growth method employed in Ref. [29|. In our 
high-quality samples, besides observing zero-resistivity, 
we were able to measure the specific-heat anomaly in the 
superconducting phase, which indicated that the super- 
conductivity in this material is a bulk feature and, more- 
over, appears to have a full energy gap^i which is a pre- 
requisite to topological superconductivity^ In this pa- 
per, we describe in detail the newly applied preparation 
method of Cu ;E Bi2Se3 to employ electrochemical interca- 
lation, and present the electronic phase diagram of T c vs 
Cu concentration x. We also show how the superconduct- 
ing shielding fraction of the electrochemically-prepared 
samples changes with x, which bears an intriguing impli- 
cation of intrinsically inhomogeneous superconductivity. 



II. SAMPLE PREPARATIONS 

Pristine Bi2Se3 has a layered crystal structure (R3m, 
space group 166) consisting of stacked Se-Bi-Se-Bi-Se 
quintuple layers. The rhombohedral [111] direction is de- 
noted as the c axis and the (111) plane as the ab plane. 
The neighboring quintuple layers are only weakly van- 
der-Waals bonded to each other. Cu is known to enter 
Bi2Se3 either as an intercalant in the van-der-Waals gaps 
or as a substitutional defect to replace Bi. In the former 
case Cu + is formed and acts as a donor, whereas in the 
latter case it creates two holes by replacing three Bi 6p 
electrons by one Cu As electron upon forming a a bond. 
Therefore Cu acts as an ambipolar dopant for Bi2Se3 J 32 i 33 
The melt-growth method employed by Hor et al. 29 can- 
not avoid the formation of substitutional Cu defects; in 
contrast, our electrochemical method has a distinct ad- 
vantage to promote only the intercalation of Cu. 

For the present experiments, we first grew single crys- 
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FIG. 1. (color online) (a) Sketch of our simplified electro- 
chemical intercalation setup. The counter electrode (CE) is 
also used as reference electrode (RE). The current direction 
is indicated, (b) Typical time dependence of the voltage drop 
measured between working and counter electrode. The inset 
summarizes the development of the diamagnetism (plotted in 
terms of the shielding fraction) in a sample with x = 0.31 
upon annealing at different temperatures. 



is a destructive method but is good at giving absolute 
numbers, to confirm that the x values determined from 
the weight change is indeed reliable, as will be described 
later. 

For a consistency check of our electrochemical interca- 
lation process, we have calculated the total charge trans- 
ferred during the intercalation. The experimentally de- 
termined Faraday number for our samples was usually 
up to 15% smaller than the theoretical Faraday constant, 
96485 C/mol. In this regard, it was reporte d 32 ! 33 that a 
simple physical attachment of Cu metal to Bi2Se3 crystal 
results in a perceptible diffusion of Cu + into the van-der- 
Waals gaps of Bi2Se3. This reaction can be described as 
the following oxidation and reduction formulas in an elec- 
trochemical viewpoint: 



xCu — >• xCu + + xe 



(i) 



tals of pristine Bi2Se3 by melting stoichiometric amounts 
of high-purity elemental shots of Bi (99.9999%) and Se 
(99.999%) at 850°C for 96 h in sealed evacuated quartz 
glass tubes, followed by a slow cooling to 550° C over 72 
h and annealing at that temperature for 48 h. Before 
the electrochemical intercalation procedure, those melt- 
grown Bi2Se3 single crystals were cleaved and cut into 
rectangular pieces. They were wound by 50 /im thick, 
bare Cu wire and acted as the working electrode (WE). 
A 0.5-mm thick Cu stick was used both as the counter 
(CE) and reference electrode (RE) in our simplified setup 
sketched in Fig. [3Ja). 

Since Hor et al£& reported that Cu x Bi2Se3 samples are 
sensitive to air, the pristine Bi2Se3 samples were trans- 
ferred into a glove box and the electrochemical intercala- 
tion of Cu was done in an inert atmosphere. To provide 
the electrical wiring required for electrochemical interca- 
lation processes as sketched in Fig. [TJa), our glove box 
was specially modified for air-tight electrical connections. 
For the intercalation we used a saturated solution of Cul 
powder (99.99%) in acetonitrile CH 3 CN. A current of 10 
jiA was applied to the electrodes for a suitable time pe- 
riod to give the desired Cu concentration x. The typical 
time dependence of the voltage drop between the stan- 
dard and working electrodes observed for our samples 
is shown in Fig. QJb). Its absolute value was usually 
between —10 and —30 mV. The whole process was con- 
trolled by a computer which was also used to determine 
the total transferred charge. 

After the intercalation had finished, the Cu concen- 
tration x in the sample was determined from its weight 
change: we measured the weight before and after the in- 
tercalation process with a high-precision balance with a 
resolution of 0.1 /ig, and the possible error in x was less 
than ±0.01. Note that this direct measurement of the 
Cu concentration is expected to be more accurate than 
most of the chemical analysis methods^ Nevertheless, 
we have also employed the inductively-coupled plasma 
atomic-emission spectroscopy (ICP-AES) analysis, which 



Bi 2 Se 3 + xe- -> (BiaSea)*". (2) 

The occurrence of electron transfer from Cu to Bi2Se3 
implies that the redox potential of Cu + /° is more neg- 
ative than that of (Bi2Se3)°/ :c_ . Since we are using a 
Cu wire to hold the sample and to make the electrical 
contact, it is possible that this wire became the source of 
an additional Cu + intercalation to the sample (s). This 
reaction should not be included into the charge balance, 
and hence the difference between the estimated and the 
theoretical Faraday number is naturally expected. This 
inference is supported by our observation that the sample 
mass increased even without applying a current, though 
at a much lower reaction rate. 

It is most likely that the electrochemical intercalation 
is induced by reduction of Bi2Se3 without reduction of 
Cu + , yielding (Cu + ) x (Bi2Se3) :c_ . This implies that dop- 
ing of electrons with the nominal fraction of x should 
take place upon Cu + intercalation. The reaction at the 
counter electrode seems to be oxidation of Cu metal 
rather than that of I - to I2 or I3", because the former 
redox potential is more negative than the latter one in 
an aqueous medium. 

It is important to note that the as-intercalated sam- 
ples do not superconduct yet. It turned out that the 
samples have to be annealed in order to establish super- 
conductivity. For this purpose, the samples were sealed 
under vacuum in quartz glass tubes and put into a muffle 
furnace. For the determination of the optimum anneal- 
ing temperature, one sample with x = 0.31 was annealed 
subsequently for 2 h at different temperatures starting 
at 500° C. After each annealing run, the diamagnetic re- 
sponse was checked, and the result is shown in the in- 
set of Fig. DJb). A trace of diamagnetism indicating the 
appearance of a superconducting phase was detected af- 
ter annealing at 530° C, and the annealing temperature 
of 560° C was found to yield the largest shielding frac- 
tion. When the sample was annealed at higher temper- 
atures than 560° C, the shielding fraction was found to 
be reduced irreversibly. (The onset T c was essentially 
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FIG. 2. (color online) EPMA analyses of the cleaved sur- 
face of a sample with x = 0.31. Panels (a) and (b) summarize 
scans along two lines of 10 /im and 1 mm lengths, respectively, 
sketched in panel (c). One can see that the local Cu concen- 
tration is uniform on the /im scale, but it is inhomogeneous 
on the sub-mm scale. 



independent of the annealing temperature.) Therefore, 
the samples used for determining the phase diagram pre- 
sented here were treated in the following way: They were 
heated up to 540° C in 1 h, and then the temperature 
was gradually increased to 560° C in 40 min to avoid any 
overheating; the samples were kept at 560° C for 2 h and 
eventually quenched by dropping the quartz glass tubes 
into cold water ^ We will come back later to the ques- 
tion of what is happening during the annealing process 
to activate the superconductivity. 

To confirm the accuracy of the x values determined 
by the mass change (and also to make sure that the Cu 
content does not change appreciably during the annealing 
process), we measured the x values of eight of our samples 
in the post-annealed state by using the ICP-AES analysis 
(the sample mass was between 15 to 38 mg). For this 
destructive analysis, the whole sample was dissolved in 
nitric acid HNO3. The Cu concentrations obtained from 
the ICP-AES analyses for the eight samples agreed with 
those obtained from the mass change within ±0.014^& 
giving confidence in the x values reported in this paper. 

In addition to the ICP-AES analysis, we have employed 
the electron-probe micro-analyzer (EPMA) to check the 




FIG. 3. (color online) (a) Temperature-dependent resistivity 
of CuxBi 2 Se 3 for x = 0.14, 0.29, and 0.40. The inset gives 
an enlargement around the superconducting transition, (b) 
Temperature dependence of the Hall coefficient Ru and the 
x dependence of the charge- carrier concentration n (inset), 
(c) Temperature-dependence of the superconducting shield- 
ing fraction in B = 0.2 mT after ZFC for various Cu concen- 
trations 0.12 < x < 0.60. For x = 0.29 the ZFC and FC data 
are exemplarily shown in the inset. 



distribution of Cu within the sample after the annealing. 
For this analysis, one sample with x = 0.31 was cleaved 
and subsequently scanned over distances of 10 fim and 1 
mm, as sketched in Fig. [2fc). The intensity of the char- 
acteristic X-ray of each element is plotted vs the scan 
position for the two different scan lengths in Figs. 2(a) 
and 2(b). One can see that the distributions of Bi and Se 
are essentially homogeneous on any length scale, as evi- 
denced by the constant intensity of the respective charac- 
teristic X-ray. On the other hand, the distribution of Cu 
shows a variation of up to ~20% on the sub-mm length 
scale [Fig. 2(b)], although on the 10- jam length scale the 
Cu distribution is usually homogeneous [Fig. 2(a)]. Since 
the averaged Cu concentration of this sample determined 
from the mass change was 0.31, the spatial variation of 
~20% corresponds to the variation in x of ~0.06. 



III. SAMPLE CHARACTERIZATIONS 

The superconducting samples were characterized by 
measuring the dc magnetization M and transport prop- 
erties. To take the magnetization data, a commercial 
SQUID magnetometer (Quantum Design, MPMS) was 
used with the magnetic field applied parallel to the ab 
plane. The samples were cooled down in zero-magnetic 
field (zero-field cooled, ZFC) to the lowest accessible tern- 



perature T = 1.8 K; then a small dc field of B = 0.2 mT 
was applied and the magnetization was measured upon 
increasing temperature. After passing through T c , de- 
fined as the onset of the drop in the M(T) curves, data 
were again taken upon decreasing the temperature back 
to 1.8 K (field cooled, FC). The superconducting shield- 
ing fraction of a sample was estimated from its magnetic 
moment at T = 1.8 K after ZFC. The resistivity p xx and 
the Hall coefficient Rh were measured by a standard six- 
probe technique, where the electrical current was applied 
in the ab plane. 

With our electrochemical intercalation technique, we 
have successfully synthesized samples which supercon- 
duct above 1.8 K for Cu concentrations of 0.09 < x < 
0.64. Figure El^a) shows the resistivity data of Cu^E^Sea 
for three selected Cu concentrations x = 0.14, 0.29, and 
0.40 measured in zero field. All those samples exhibit a 
metallic temperature dependence above T c and show zero 
resistance below T c , see the expanded view in the inset of 
Fig. [3fa). With increasing Cu concentration, the abso- 
lute value of p xx increases; especially, between x = 0.29 
and 0.40 a strong rise in the absolute value is observed, 
which implies that a high Cu concentration enhances the 
disorder in the samples. FigureJ^b) summarizes the tem- 
perature dependences of Rh for the three samples, which 
are generally weak. The inset shows the ^-dependence of 
the charge-carrier concentration n determined from the 
low-temperature value of Rh- It is striking that, despite 
the factor of three difference in the Cu concentration, the 
change in n is very small: n = 1.5 xlO 20 , 1.3 xlO 20 , and 
1.7xl0 20 cm- 3 for x = 0.14, 0.29, and 0.40, respectively. 
Furthermore, those values correspond to only ~2% of 
electron doping. This is totally inconsistent with the ex- 
pectation that there should be electron carriers with the 
nominal fraction of x in Cu^E^Ses. Therefore, there 
must be some side reaction taking place in Cu^E^Ses to 
significantly reduce the actual electron carriers. In this 
regard, a similar problem was previously noted in a study 
of Cu intercalation into Bi2Te3, and it was proposed that 
Cu + reacts with the matrix to form a four-layer lamellar 
structure Cu-Te-Bi-Te that annihilates electrons^ If a 
similar reaction occurs in Cu^E^Ses, it would be 



Cu + + 2Bi 2 Se 3 



[CuBiSe 2 ][Bi 3 SeJ] + /i. (3) 



The formation of this type of structural defects can in- 
deed explain the small n values observed in Cua3i2Se3. 

Figure [3](c) summarizes the temperature dependences 
of the shielding fraction in samples with various x val- 
ues. The onset temperatures of superconductivity in the 
magnetization data are slightly lower than the onset of 
the resistivity transition in respective samples, which is 
usual for disordered superconductors. The right inset 
of Fig. E^c) shows the ZFC and FC data for x = 0.29, 
and other samples exhibited similar differences between 
the two. We note that samples with x > 0.70 were also 
prepared, but they did not show any superconducting 
transition above 1.8 K; however, it is possible that they 
exhibit superconductivity at lower temperatures. 
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FIG. 4. (color online) Cu concentration x dependence of (a) 
the critical temperature and (b) the shielding fraction at 1.8 
K. These T c data were extracted from the onset of super- 
conductivity in magnetization measurements, the shielding- 
fraction data from the ZFC magnetization at T — 1.8 K, see 
Fig. He). 



IV. PHASE DIAGRAM AND DISCUSSIONS 

Figure Ufa) shows the phase diagram of T c vs x based 
on our measurements of more than 40 samples, and 
Fig. [U(b) summarizes the shielding fraction of those sam- 
ples at 1.8 K. The lowest x value at which we found 
superconductivity was 0.09, where T c was already 3.45 
K; however, in spite of this relatively high T c , the super- 
conducting signal was very weak with a shielding fraction 
of only around 1%. For slightly larger x values around 
0.12-0.15 the maximum values of T c are found, but the 
shielding fractions were less than 20%, which is in agree- 
ment with the earlier report by Hor et alS^ Intriguingly, 
T c gradually decreases with increasing x down to T c w 2.2 
K for x = 0.64, which is the highest Cu concentration for 
which we could clearly confirm the superconductivity. 

Although T c seems to be quite robust and reproducible, 
the shielding fraction exhibits a scattering among differ- 
ent samples even with similar x values, as can be seen in 
Fig. Hfb). Nevertheless, there is a clear trend that sam- 
ples with maximal shielding fractions are obtained only 
for 0.3 < x < 0.5 and the shielding fraction is reduced 
systematically (almost linearly) as x becomes smaller 
than x ~ 0.3. This trend, together with the EPMA re- 
sult that the actual Cu concentration varies on a sub-mm 
scale, suggests that Cu x Bi2Se3 has a tendency to phase- 
segregate and that there is a spontaneous formation of 
small islands in which the local Cu concentration is essen- 
tially the optimum value for superconductivity. In view 
of the fact that the T c of this material is uncharacteristi- 
cally high for a low-carrier-density superconductor 38 with 
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n ~ 10 20 cm -3 , the present observation seems to imply 
that the inhomogeneous nature of the superconductivity 
may partly be responsible for the anomalously high T c . 39 
Furthermore, the unexpected drop of the shielding frac- 
tion for x > 0.5 suggests that the inhomogeneity may be 
necessary for the occurrence of superconductivity in this 
material. 



If Cu^E^Sea indeed has a tendency toward phase seg- 
regation to spontaneously form small superconducting is- 
lands, it naturally explains why a high- temperature an- 
nealing is necessary for activating the superconductivity: 
It allows Cu atoms to move and promotes the forma- 
tion of superconducting islands. If annealed at too high 
a temperature, the material seems to melt partly, the 
crystal structure degrades and eventually the supercon- 
ductivity is destroyed irreversibly. Indeed, the melting 
temperature of Cu x Bi2Se3 appears to be systematically 
reduced with x, and our test sample with x = 1.4, which 
was initially rectangular shaped, became sphere-like af- 
ter annealed at 560° C. (Note that the melting point of 
pure Bi2Se3 is 706°C). Therefore, the optimal annealing 
temperature corresponds to the one that maximally pro- 
motes the motion of Cu atoms while avoiding degradation 
of the matrix due to partial melting. 



V. SUMMARY 

In summary, we have synthesized superconducting 
Cu^E^Sea by an electrochemical intercalation method, 
which yields samples with much higher superconducting 
shielding fractions compared to those prepared by the 
melt-growth technique previously employed. The super- 
conductivity was observed for 0.09 < x < 0.64, which 
is much wider than previously reported, and the transi- 
tion temperature T c was found to show an unexpected 
monotonic decrease with x. Also, the largest attainable 
shielding fraction was found to be strongly dependent on 

and its x dependence suggests that Cu^E^Sea has a 
tendency to phase-segregate and form small supercon- 
ducting islands. 
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minimize the sample preparation time. 
The detailed results of the ICP-AES analyses are the fol- 
lowing: the x values of the eight samples determined from 
the mass change were 0.117, 0.128, 0.266, 0.286, 0.293, 
0.295, 0.308, and 0.318; for these samples, the x values 
determined by the ICP-AES analyses were 0.114, 0.128, 
0.279, 0.300, 0.299, 0.282, 0.303, and 0.319, respectively; 
therefore, the differences between the two were 0.003, 

0. 000, -0.013, -0.014, -0.006, 0.013, 0.005, and -0.001. 
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